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Measurements of ip and tp' production from experiment NA50 at the CERN SPS are 
compared to calculations based on a hadronic model of charmonium suppression devel- 
oped previously. Data on centrality dependence and total cross sections are in good 
accord with these predictions. Uncertainties in theoretical quantities such as NA50's L 
parameter are discussed. 

1. Introduction 

Has the quark gluon plasma been discovered at the CERN SPS? Experiment NA50 
has reported an abrupt decrease in -0 production in Pb+Pb collisions at 158 GeV 
per nucleon p]]. Specifically, the collaboration presented a striking 'threshold effect' 
in the ^-to-continuum ratio by plotting it as a function of a calculated quantity, 
the mean path length of the through the nuclear medium, L, as shown in fig. la. 
This apparent threshold has sparked considerable excitement as it may signal de- 
confinement in the heavy Pb+Pb system . 

In this talk I report on work with Ramona Vogt in ref. Q comparing Pb results 
to predictions |^] using a hadronic model of charmonium suppression. We first 
demonstrate that the behavior in the NA50 plot, fig. la, is not a threshold effect but, 
rather, reflects the approach to the geometrical limit of L as the collisions become 
increasingly central. When plotted as a function of the measured neutral transverse 
energy Et as in fig. lb, the data varies smoothly as in S+U measurements in fig. 3b 
below [Q, ||, 0, ||, [|. The difference between S+U and Pb+Pb data lies strictly 
in the relative magnitude. To assess this magnitude, we compare i\) and tp' data 
to expectations based on the hadronic comover model 0, pi. The curves in fig. 1 
represent our calculations using parameters fixed earlier in Ref. Our result is 
essentially the same as the Pb+Pb prediction in Q . 

Our primary intention is to demonstrate that there is no evidence for a strong 
discontinuity between pA, S+U and Pb+Pb data. However, to quote Maurice 
Goldhaber, ". . . absence of evidence is not evidence of absence." Our secondary 
goal is to show that our model predictions agree with the new Pb+Pb data. The 
consistency of these predictions is evident from the agreement of our old pA and S+U 
calculations with more recent NA38 and NA51 data. Nevertheless, the significance 
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Figure 1: (a) The NA50 § comparison of ip production in Pb+Pb and S+U colli- 
sions as a function of the average path length L, see eq. (3). B is the ip — * /i + [f~ 
branching ratio, (b) Transverse energy dependence of Pb+Pb data. Curves in (a) 
and (b) are computed using eqs. (4-6). 



of this result must be weighted by the fact that all pA and AB data are preliminary 
and at different beam energies. 

In this work, we do not attempt to show that our comover interpretation of the 
data is unambiguous - this is certainly impossible at present. 

2. Nucleons and Comovers 

The hadronic contribution to charmonium suppression arises from scattering 
of the nascent tp with produced particles - the comovers - and nucleons ||, 
To determine the suppression from nucleon absorption of the ip, we calculate the 
probability that a cc pair produced at a point (b, z) in a nucleus survives scattering 
with nucleons to form a ip. The standard |}| [l0| result is 

dz p A (b, z)cr^ N } (1) 

where pa is the nuclear density, b the impact parameter and ct^n the absorption 
cross section for ^-nucleon interactions. One can estimate Sa ~ exp{— u^n p$La\, 
where La is the path length traversed by the cc pair. 

Suppression can also be caused by scattering with mesons that happen to travel 
along with the cc pair (see refs. in Q). The density of such comovers scales roughly 
as Et- The corresponding survival probability is 
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where n is the comover density and t is the time in the ip rest frame. We write 
S co ~ exp{ — [3Et}, where [3 depends on the scattering frequency, the formation 
time of the comovers and the transverse size of the central region, Rt, cf. eq. (8). 

To understand the saturation of the Pb data with L in fig. la, we apply the 
schematic approximation of Rcf. p0[ for the moment to write 



< + B u-(Et) 



t \'\ oc (SaSbSco) ~ e -^P" L e-^, (3) 



where the brackets imply an average over the collision geometry for fixed Et and 
<t{Et) = da/dET- The path length L = (La + Lb) and transverse size Rt de- 
pend on the collision geometry. The path length grows with Et, asymptotically 
approaching the geometric limit Ra + Rb- Explicit calculations show that nucleon 
absorption begins to saturate for b < Ra, where Ra is the smaller of the two nuclei, 
see fig. 4 below. On the other hand, Et continues to grow for b < Ra due, e.g., 
to fluctuations in the number of NN collisions. Equation (2) falls exponentially in 
this regime because (3, like L, saturates. 

In fig. lb, we compare the Pb data to calculations of the ?/>-to-continuum ratio 
that incorporate nucleon and comover scattering. The contribution due to nucleon 
absorption indeed levels off for small values of b, as expected from eq. (3). Comover 
scattering accounts for the remaining suppression. 

These results are predictions obtained using the computer code of Ref. with 
parameters determined in Ref. ||. However, to confront the present NA50 analysis 
[Q, we account for changes in the experimental coverege as follows: 

• Calculate the continuum dimuon yield in the new mass range 2.9 < M < 
4.5 GeV. 



Adjust the E T scale to the pseudorapidity acceptance of the NA50 calorimeter, 
l.K r) < 2.3. 



The agreement in fig. 1 depends on these updates. 
3. J/tjj Suppression 

We now review the details of our calculations, highlighting the adjustments as 
we go. For collisions at a fixed b, the V'-production cross section is 

O) = <J$ N Jd 2 sdzdz' PA (s, z) PB ib - s, z') S, (4) 

where S = SaSbS co is the product of the survival probabilities in the projectile A, 
target B and comover matter. The continuum cross section is 

= o-^Jd 2 sdzdz'p A (s,z)p B (b - s,z'). (5) 

The magnitude of (4,5) and their ratio are fixed by the elementary cross sections 
a^ N and <T ^+^- • We calculate a^ N using the phenomenologically-successful color 
evaporation model p[. The continuum in the mass range used by NA50, 2.9 < 
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M < 4.5 GeV, is described by the Drell-Yan process. To confront NA50 and NA38 
data in the appropriate kinematic regime, we compute these cross sections at leading 
order following [y|, |l| using GRV LO parton distributions with a charm i^-factor 
K c = 2.7 and a color evaporation coefficient F.,p = 2.54% and a Drell-Yan if-factor 
Kjyy = 2.4. Observe that these choices were fixed by fitting pp data at all available 
energies pTT[ | . Computing f° r 2-9 < M < 4.5 GeV corresponds to the first 

update. 

To obtain Et dependent cross sections from eqs. (4) and (5), we write 

a AB (E T ) = J d 2 bP(E T ,b)a AB (b). (6) 

The probability P(Et, b) that a collision at impact parameter b produces transverse 
energy Et is related to the minimum-bias distribution by 

Cm in (Et) = J d 2 bP(E T ,b). (7) 

We parametrize P(E T ,b) = Cexp{-(E T -E T ) 2 /2A}, where E T (b) = ej\f(b), 
A(b) = ueET(b), C(b) = (27rA(fe)) _1 and Af(b) is the number of participants (see, 
e.g., Ref. |Q). We take e and u> to be phenomenological calorimeter-dependent 
constants. 

We compare the minimum bias distributions for total hadronic Et calculated 
using eq. (7) for e = 1.3 GeV and w = 2.0 to NA35 S+S and NA49 Pb+Pb data 
[ O) . The agreement in fig. 2a builds our confidence that eq. (7) applies to the 
heavy Pb+Pb system. Figure 2b shows the distribution of neutral transverse energy 
calculated using eqs. (5) and (6) to simulate the NA50 dimuon trigger. We take 
e = 0.35 GeV, co — 3.2, and a N + _ ~ 37.2 pb as appropriate for the dimuon- 
mass range 2.9 < M < 4.5 GeV. The E T distribution for S+U -> fi + /.r + X 
from NA38 was described || using e = 0.64 GeV and w = 3.2 - the change in 
e corresponds roughly to the shift in particle production when the pseudorapidity 
coverage is changed from 1.7 < rj < 4.1 (NA38) to 1.1 < 77 < 2.3 (NA50). Taking 
e = 0.35 GeV for the NA50 acceptance is the second update listed earlier. We 
now apply eqs. (1,2,4) and (5) to charmonium suppression in Pb+Pb collisions. To 
determine nucleon absorption, we used pA data to fix ct^n ~ 4.8 mb in Ref. ||. 
This choice is in accord with the latest NA38 and NA51 pA data, see fig. 3a. To 
specify comover scattering [^), we assumed that the dominant contribution to if) 
dissociation comes from exothermic hadronic reactions such as p + if) — > D + D. We 
further took the comovers to evolve from a formation time tq ~ 2 fm to a freezeout 
time tf ~ Rt/v 1c i following Bjorken scaling, where v Te i ~ 0.6 is roughly the average 
if) — p relative velocity. The survival probability, eq. (2), is then 

= exp{-a co v Te in T \n(R T / v ic1 t q )} (8) 

where a co ~ 2ct^at/3, Rt ~ Ra and hq is the initial density of sufficiently massive 
p, u) and 77 mesons. To account for the variation of density with Et, we take no = 
hqEt/ Et(0) 0- A value no — 0.8 fm~ 3 was chosen to fit the central S+U datum. 
Since we fix the density in central collisions, this simple ansatz for S co may be 
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Figure 2: Transverse energy distributions from eq. (7). The S-Pb comparison (a) 
employs the same parameters. 



inaccurate for peripheral collisions. [Densities ~ 1 fm typically arise in hadronic 



models of ion collisions, e.g., refs. 15 . The internal consistency of hadronic models 
at such densities demands further study] 

We expect the comover contribution to the suppression to increase in Pb+Pb 
relative to S+U for central collisions because both the initial density and lifetime 
of the system can increase. To be conservative, we assumed that Pb and S beams 
achieve the same mean initial density. Even so, the lifetime of the system essentially 
doubles in Pb+Pb because Rt ~ Ra increases to 6.6 fm from 3.6 fm in S+U. 
The increase in the comover contribution evident in comparing figs, lb and 3b 
is described by the seemingly innocuous logarithm in eq. (8), which increases by 
w 60% in the larger Pb system. 

In Rcf. [El, we pointed out that comovers were necessary to explain S+U data 
from the NA38 1991 run 0. Data just released Q| from their 1992 run support this 
conclusion. The '91 ip data were presented as a ratio to the dimuon continuum in the 
low mass range 1.7 < M < 2.7 GeV, where charm decays are an important source of 
dileptons. On the other hand, the '92 ip data jl], [| are given as ratios to the Drell- 
Yan cross section in the range 1.5 < M < 5.0 GeV. That cross section is extracted 
from the continuum by fixing the X -factor in the high mass region [ |l4[ . To compare 
our result from Ref. j^j to these data, we scale the '92 data by an empirical factor. 
This factor is » 10% larger than our calculated factor (92) /a™ t (91) » 0.4; 
these values agree within the NA38 systematic errors. [NA50 similarly scaled the 
'92 data to the high-mass continuum to produce fig. la.] Because our fit is driven 
by the highest Et datum, we see from fig. 3b that a fit to the '92 data would not 
appreciably change our result. Note that a uniform decrease of the ratio would 
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Figure 3: (a) pA cross sections |lj in the NA50 acceptance and (b) S+U ratios from 
'91 @ and '92 runs. The '92 data are scaled to the '91 continuum. The dashed 
line indicates the suppression from nucleons alone. The pp cross section in (a) is 
constrained by the global fit to pp data in ref. pT[ . 



increase the comover contribution needed to explain S+U collisions. 

NA50 and NA38 have also measured the total ■i/'-production cross section in 
Pb+Pb |i] and S+U reactions To compare to that data, we integrate eqs. (4, 6) 
to obtain the total (a/AB)^ = 0.95 nb in S+U at 200 GeV and 0.54 nb for Pb+Pb 
at 158 GeV in the NA50 spectrometer acceptance, 0.4 > xf > and —0.5 < cos 9 < 
0.5 (to correct to the full angular range and 1 > xf > 0, multiply these cross sections 
by » 2.07). The experimental results in this range are 1.03 + 0.04 + 0.10 nb for S+U 
collisions [§ and 0.44 ± 0.005 ± 0.032 nb for Pb+Pb reactions J|. Interestingly, in 
the Pb system we find a Drell-Yan cross section (a/AB) DY = 37.2 pb while NA50 
finds (a/AB) DY = 32.8 ± 0.9 ± 2.3 pb. Both the ip and Drell-Yan cross sections 
in Pb+Pb collisions are somewhat above the data, suggesting that the calculated 
rates at the NN level may be ~ 20 — 30% too large at 158 GeV. This discrepancy is 
within ambiguities in current pp data near that low energy |pd| . Moreover, nuclear 
effects on the parton densities omitted in eqs. (4,5) can affect the total S and Pb 
cross sections at this level. 

We remark that if one were to neglect comovers and take ct^jv = 6.2 mb, one 
would find (a/AB)^ = 1.03 nb in S+U at 200 GeV and 0.62 nb for Pb+Pb at 
158 GeV. The agreement with S+U data is possible because comovers only con- 
tribute to the total cross section at the ~ 18% level in the light system. This is 
expected, since the impact-parameter integrated cross section is dominated by large 
b and the distinction between central and peripheral interactions is more striking 
for the asymmetric S+U system. As in Ref. M, the need for comovers is evident 



6 



15 



10 - 



a 



5 - 







1 1 I 1 1 I 

; L(E T ) from [1] 


1 1 1 1 1 1 1 1 1 1 

L 


L 




















^ <b> 


\< b > 


■ S+U " ^ 


- Pb+Pb v ■ 



20 



50 

NA38 



80 50 



150 



100 

■"p _ . , ^ T NA5 ° (GeV) 
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impact parameter (b) (dot-dashed). The solid line covers the measured Et range. 
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for the -ET - dependent ratios, where central collisions are singled out. 
4. Saturation and the Definition of L 

To see why saturation occurs in Pb+Pb collisions but not in S+U, we compare 
the NA50 L(Et) m to the average impact parameter (t»}(£+) in fig. 4. To best un- 
derstand fig. la, we show the values of L(Et) computed by NA50 for this figure. We 
use our model to compute (b) = (bT AB ) / (T AB ) , where (/(&)) = Jd 2 b P(E T ,b)f(b) 
and T AB = JcPsdzdz' p A (s, z)p B (b — s,z'). [Note that NA50 reports similar val- 
ues of (&)(£'t) JH-] In the Et range covered by the S experiments, we see that 
(6) is near ~ R$ = 3.6 fm or larger. In this range, increasing b dramatically re- 
duces the collision volume and, consequently, L. In contrast, in Pb+Pb collisions 
(b) <C i?pb = 6.6 fm for all but the lowest Et bin, so that L does not vary appre- 
ciably. 

To understand the sensitivity of fig. la to the definition of the path length, we 
now estimate L(Et) [lq| - We identify (3) with the exact expression formed from 
the ratio of (4) and (5). Expanding in <t^n and neglecting comovers, we find: 

L(E T ) = {2p (T AB )}- 1 l^J d 2 s [T A {s)] 2 T B {b -s) + [T B (b - s)] 2 T A ( S )^ , (9) 

where T A (s) — J p A (s,z)dz. In fig. 5 we compare the NA50 L(Et) to the path 
length calculated using two assumptions for the nuclear density profile: our re- 
alistic three-parameter Fermi distribution and the sharp-surface approximation 
p = po<d(R A — r). NA38 Jl7| obtained L for S+U using the empirical prescription 
of ref . Jl(| , while NA50 calculated L assuming the sharp-surface approximation JlJ] . 
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Figure 5: NA50 L{Et) [1] (points) compared to calculations for realistic nuclear 
densities (solid), as used here, and for a sharp-surface approximation (dot-dashed). 
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Figure 6: NA50 data replotted with a realistic L{Et) from (9). 
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[0 | and (b) NA38 S+U data §] (filled points) and preliminary data @. 



Indeed, we see that the NA50 Pb+Pb values agree with our sharp-surface result, 
while the NA38 S+U values are nearer to the realistic distribution. 

To see how the value of the path length can affect the appearance of fig. la, 
we replot in fig. 6 the NA50 data using L(Ej<) from (9) with the realistic density. 
We learn that the appearance of fig. la is very sensitive to the definition of L. 
Furthermore, with a realistic L, one no longer gets the impression given by the 
NA50 figure |J of Pb+Pb data "departing from a universal curve." Nevertheless, 
the saturation phenomena evident in fig. la does not vanish. Saturation is a real 
effect of geometry. 

5. ip' Suppression 

To apply eqs. (4-6) to calculate the ip'-to-ip ratio as a function of Er, we 
must specify <J^, N , ct^'at, and cr^> co . Following Ref. fUjl , we use pp data to fix 
Ba$, N /B<7$ N = 0.02 (this determines ify). The value of a^i n depends on whether 
the nascent ip' is a color singlet hadron or color octet cc as it traverses the nucleus. In 
the singlet case, one expects the absorption cross sections to scale with the square of 
the charmonium radius. Taking this ansatz and assuming that the ip' forms directly 
while radiative x decays account for 40% of ip production, one expects a^i ~ 2.1(7,/, 
for interactions with either nucleons or comovers [pi. For the octet case, we take 
a^p'N ~ c-ipN and fix a^> co w 12 mb to fit the S+U data. In fig. 7a, we show that 
the singlet and octet extrapolations describe pA data equally well. 

Our predictions for Pb+Pb collisions are shown in fig. 8. In the octet model, the 
entire suppression of the ip'-to-ip ratio is due to comover interactions. In view of the 
schematic nature of our approximation to 5 co in eq. (8), we regard the agreement 
with data of singlet and octet extrapolations as equivalent. 




Figure 8: Comover suppression in Pb+Pb — > ip' + X. 



6. Summary 

In summary, the Pb data Q cannot be described by nucleon absorption alone. 
This is seen in the NA50 plot, fig. la, and confirmed by our results. The saturation 
with L but not Et suggests an additional density-dependent suppression mecha- 
nism. Earlier studies pointed out that additional suppression was already needed to 
describe the S+U results ||; recent data Q support that conclusion (see, however, 
[||). Comover scattering explains the additional suppression. Nevertheless, it is un- 
likely that this explanation is unique. SPS inverse-kinematics experiments (B < A) 
and AGS pA studies near the ip threshold can help pin down model uncertainties. 

After the completion of || , several cascade calculations |l5| have essentially con- 
firmed our conclusions. This confirmation is important, because such calculations 
do not employ the simplifications (e.g. no oc Et) needed to derive (8). In par- 
ticular, these models calculate Et and the comover density consistently. Some of 
these authors took ct^jv ~ 6 mb (instead of ~ 5 mb) to fit the NA51 data in fig. 3a 
somewhat better. 

I am grateful to Ramona Vogt for her collaboration in this work. I also thank 
C. Gerschel and M. Gonin for discussions of the NA50 data, and M. Gyulassy, 
R. Pisarski and M. Tytgat for insightful comments. 
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